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ABSTRACT 

Classical Cepheids are primary distance indicators playing a fundamental role 
in the calibration of the extragalactic distance scale. The possible dependence 
of their characteristic Period-Luminosity (PL) relation on chemical composition 
is still debated in the literature, and the behaviour of these pulsators at very 
low metallicity regimes is almost unexplored. In order to derive constraints on 
the application of the Period-Luminosity relation at low metal abundances, we 
investigate the properties of the few ultra-low metallicity (Z « 0.0004) Cepheids 
recently discovered in the Blue Compact Dwarf galaxy IZwl8. To this purpose we 
have computed an updated and extended set of nonlinear convective models for 
Z = 0.0004 and Y = 0.24, spanning a wide range of stellar masses, and taking 
into account the evolutionary constraints for selected luminosity levels. As a 
result we are able to predict the topology of the instability strip, the variations 
of all the relevant quantities along the pulsation cycle, including the morphology 
of the light curves, the theoretical Period-Luminosity-Color, Period- Wesenheit 
and Period-Luminosity relations at such a low metallicity. For each of these 
relations we provide the appropriate coefficients for fundamental mode pulsators 
with Z=0.0004. By comparing these results with the properties of more metal 
rich Cepheids we find that the synthetic PL relations for Z = 0.0004 are steeper 
than at higher Z, but similar to the Z = 0.004 ones, thus suggesting a leveling 
off of the metallicity effect towards the lowest Zs. 

Subject headings: galaxies: dwarf — galaxies: irregular — galaxies: individual 
(I Zw 18) — galaxies: stellar content — galaxies: Cepheids — Distance scale 
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1. Introduction 

With their characteristic Period-Luminosity (PL) and Period-Luminosity-Color (PLC) 
relations, Classical Cepheids play a relevant role as distance indicators for the definition of 
the extragalactic distance scale. The application of a Large Magellanic Cloud (LMC)-based 
PL relation to external galaxies observed with the Hubble Space Telescope (HST), has led 
to the calibration of s econdary distance indicators and in tur n to an estimate of the Hubble 



constant (H , see e.g. iFreedman et al.ll200ll ; ISaha et al.l 120011 ) . 



The theoretical explanation for the observational evidence of a PL relation for Classical 
Cepheids relies on the assumption that intermediate mass stars undergoing central Helium 
burning are characterized by a Mass-Luminosity (ML) relation, as predicted by stellar 
evolution models. 

The ML relation is significantly dependent on several physical and numerical 
assumptions adopted in stellar evolution models, thus the comparison between observations 
and Cepheid evolutionary/pulsation models provides a unique insight into the stellar 
evolution and pulsation physics. 

An important issue still under debate is the universality of the Cepheid PL relations, 
and, in turn, the possibility of applying LMC-calibrated Cepheid PL relations to infer the 
distance of any galaxy containing Cepheids, independently of its chemical composition. In 
the last decade many efforts were devoted to investigate the dependence of the Cepheid 
properties on metallicity, since a metallicity effect could produce significant systematic 
errors in the evaluation of the extragalactic distance scale and in turn on Hq. However, no 
general consensus has been reached yet. 

On the theoretical side, linear non adiabatic models mostly suggest that a 
variation in chemical composition produces negligible effects on the PL relations 
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998 



Sandage. Bell, fc Tripicco 



Bon o. Marconi, fc Stellingweri 



see e .g. IChiosi. Wood, fe Capitaniol Il993t lAlibert et al.l Il999l ; ISaio &: Gautschy 



1999 



1999) , while nonlinear convective pulsation models 



Fiorentino et al 



2002 



Marconi. Musella. fc Fiorentino 



20051 ) predict a significant metallicity effect on the Cepheid instability strip (here i nafter 



IS) topology and on t 



he P 



and 



Fiorentino et al. 



Caputo. Marconi, fc Musellal f l2000l ) 



relations. For instance, 
( 120021 ) find that the synthetic PL relations become shallower as the 
metallicity increases, and the size of the effect depends also on the photometric band, 
with a decreasing sensitivity as wavelength increases from optical to near infrared (NIR) 
bands. Consequently, metal-rich pulsators with periods longer than five days have fainter 
optical magnitudes than metal-poor pulsators. This scenario is further complicated by 
the theoretical finding that the helium mass fraction Y also plays a role at the highest 
me tallicities (Z > . 02), with the slope of the linear PLs increas ing as Y increases at fixed 



Z ( IFiorentino et al. 



2002 



Marconi. Musella. fc Fiorentino 



20051 ). In turn, the metallicity 



correction to the predicted distance moduli varies with the assumed AY/AZ, showing a sort 
of turnover around solar metallicity. In particular, for P > 20 days and [ O/H] > 0.2 dex. 



as m easured in several spiral galaxies observed by the HST Key Project (IFreedman et al. 



20011 ) . the average predicted metallicity correction varies from ~ —0.2 mag to ~ +0.25 



mag, as AY/AZ increases from 2.0 to 3.5. 

On the empirical side, several authors suggest that metal-rich Cephei ds are, at fixed 



period, brighter than metal-poor ones, eithe r over the entire period range ( 



1998 



Kanbu 



Macri et al. 



r et al. 



2003 



Storm et al. 



2004 



Groenewegen et al. 



2004 



Kennicutt et al 



Sakai et al 



2006) , or at least 



Sandage. Tammann. fc Reindl 



or p eri ods shorter than ~ 25 days (jSasselov et al. 



2004; 



1997 



20041 ). 



Sakai et al. 



(120041 ) by comparing distances based 



on Cepheids and on the Tip of the Red Giant Branch (TRGB) method for a selected 
sample of spiral galaxies, found that the metallicity effect, usually expressed in the form 
7 = 6fio/5\ogZ, is equal to —0.25 mag dex -1 . However, the revised TRGB distances 
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by iRizzi et al.l (120071 ) no longer support a 7 valu e of —0.25 



Bono et al 



mag dex 1 and are in 



20081 . for details). Other 



better agreement with the theoretical results (see 

em pirical results tha t seem to su pport the nonlinear theoretical scenario were obtained 
by 



Romaniello et al. 



(2005 



20081 ) . on the basis of spectroscopic [Fe/H] measurements 



for Galactic and Magellanic Cloud Cepheids. U sing only their spectrosc opic abundances 



Romaniello et al. 



(120051 ) found that the 



for Galactic Cepheids with published distances, 
metallicity correction is in qualitative agreement with the model results, and shows the 
same kind of turnover around solar metallicity as predicted by nonlinear pulsation theory. 



More recently, using direct high re solution metallici 



Galactic and Magellanic Cepheids, 



; y mea surements for a total of 68 



Romaniello et al. 



(120081 ) found that metallicity affects 



the V^-band Cepheid PL relation, with metal-rich Cepheids appearing to be sy stematically 



fainte r tha n metal-poor ones, in agreement wit h the theoretical predictions by 



(119991 ) and 



Bono et al. 



Caputo. Marconi. &: Musellal ( 120001 ). and independently of the adopted distance 



scale for Galactic and LMC Cepheids. 

On the other hand, direct empirical tests of the metallici ty effect based on the 



observatio n of Cepheids in t he outer and inner field s of M101 ( 



NGC4258 (IMacri et al. 



20061 ) . and, recently, in M33 (IScowcroft et al 



Kennicutt et al. 



1998|), in 



2009), should be taken 



with caution. Indeed, in M101 the prese nce of blended Cepheids could have affected the 



results of the test (see 



Macri et al. 



20061 ). while in the case of NGC4258 there is evidence 



against the assumed metallicity gradient, as both th e comparison wit h pulsation models 



and the most recent HII abundance measurements (IDiaz et al. 



2000) suggest a rather 



constant LMC-li ke metal-content for t he Cepheids observed i n the two selected fields (see 



Bono et al. 



by 



2008 



, for de tails). In M33, 



Scowcroft et al. 



Magrini et al 



(I2009h adopted the [O/H] gradient 



(120071 ) and conclude that blending cannot account for t 



re difference in 



distan ce modulus between the two fields. However, as also pointed out by 



Romaniello et al. 



(120081 ). these results are always based on indirect measurements of the metallicity, assumed 
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to be that corresponding to the oxygen nebular abundances of HII regions at the same 
Galactocentric distance of the Cepheids. 

Finally, there are also recent empirical results in favour of the universality of the PL 
slope: distance estim ates based on the "near-infrared surface brightn ess" (ISB) technique 



(IFouque et al. 



20071 . and references therein) and the HST parallaxes ([Benedict et al. 



20071 ) 



for ten Galactic Cepheids seem to indicate a vanishing metallicity effect between Galactic 



and Ma gellanic Cepheids. Simi 



(see e.g. 



Pietrzynski et al. 



2006 



ar suggestions have been 



2007 



Gieren et al. 



put forward by several authors 



20081 ) for Cepheids in galaxies with 



metallicities significantly lower than the LMC, but no theoretical predictions are available 
in the literature for the metallicity effect on the properties of Cepheids in this range of 
abundances or at lower Z. 

The identification and study of Cepheids in very metal-poor galaxies, is a challenging 
perspective. These galaxies represent the closest analog to primordial galaxies in the 
early universe and, as such, offer the best place where to study star formation and stellar 
evolution in a (almost) pristine environment. Some of these galaxies have also been claimed 
to be primordial systems just forming in the local Universe due to lack of detection of 
faint red stars on the red giant branch tip (TRGB, stars older than ~1 Gyr). However, 
the lack of these stars may also be due to a much larger distance of these systems that 
does not allow to resolve the faint TRGB stars. A precise determination of the distance 
of these galaxies through Cepheids allows us to verify or confute their nature of young 
nearby galaxies. In the context of the HST ACS GO program 10586 (PI: A. Aloisi) we 
have identified three Classical Cepheids in the Blue Compact Dwarf galaxy IZwl8, having 
respectively periods of 8.71, 125 and 130.3 days. A new estimate of the distance, and, in 
turn, of the age of the galaxy stellar populations were obtained both from the Classical 
Cepheids and from the galaxy's TRGB, which allowed us to rule out the possibility that 
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IZwl 8 is a truly primordial galaxy of recent formation in the local universe (jAloisi et al. 
20071 ). At the same time this pr oject also provided a fi rst sample of ultra- low metallicity 



1/5OZ ) Classical Cepheids (jFiorentino et al. 



20101 ) . thus representing an important 



benchmark for the pulsation models. 

In this paper we investigate theoretically the pulsation properties of Cepheids at this 
very low metallicity (Z=0.0004), in order to probe the metallicity effect at Z < 0.004, and 
to provide a theoretical scenario for the Classical Cepheids recently detected in IZwl8 



(jAloisi et al. 



2007 



Fiorentino et al. 



20101 ) . The organization of the paper is as follows: in 
Sect. 2 we present the new set of pulsation models; in Sections 3 and 4 we show our results 
for the topology of the IS and the behaviour of the light curves, while the multi-filter PLC, 
Period- Wesenheit and PL relations are presented and discussed in Sect. 5. Finally, in Sect. 
6 we discuss the comparison with the Cepheids observed in IZwl8 and the implications for 



the Cepheid distance scale. Conclusions are summarized in Sect. 7. 



2. Pulsation models for ultra-low metallicity Cepheids 



Using the same pulsation code (see 



Bono. Marconi, fc Stellingweri 



1999 



for details) adopted in our previous investigations of more metallic Cepheids 



Bono. Marconi. & Stellinewerf 


1999: 


Fiorentino et al. 


2002 


Marconi. Musella. & Fiorentino 



20051 ). we investigated for the first time the full amplitude behaviour of metal-poor Cepheid 
models. To this purpose we fixed the metal and helium abundances at Z = 0.0004, 
Y = 0.24 and selected model stellar masses in the range from 4 to 13 M@. Two luminosity 
levels (see Table 1) were assumed for most of the selected stellar masses corresponding to 
canonical and mild overshooting (noncanonical) evolutionary predictions, respectively. The 



canonical scenario is provided by the Pisa Evolutionary Library (PEL) for M < 11M G 
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and discussed in 



Castellani et al 



(120031 ): 



Cariulo. Degrinnocenti. fc Castellani! (]2004i )l 7 



while the noncanonical luminosity level is obtained b y adding 0.25 dex (in log L/L^) to 



Chiosi. Wood, fc Capitaniol f ll993[ ). 



the canonical value, according to the prescriptions by 
For each mass and luminosity, an extensive range of effective temperatures was explored 
for pulsation instability, in order to derive the location of the blue and red edges of the 
IS in the Hertzsprung-Russell (HR) diagram. For canonical models, we also explored 
the effect of varying the mixing length parameter l/H p , adop ted in the code to close th e 



nonlinear system of dynamica 



Bono. Marconi, fc Stellingwerl 



and convective equations (see 



1999 



Bono & Stellingwerl 



1994 



for details), in the range from the standard value of 
1.5 to 2.0. In the following two sections we discuss the results obtained for the topology of 
the IS, and for the morphological properties of the bolometric light curves, respectively. 



3. Predicted IS for Cepheids with Z=0.0004 

For each selected mass and luminosity, the effective temperatures of the hottest 
pulsating fundamental model, increased by 50 K (fundamental blue edge, hereafter FBE), 
and of the coolest pulsating fundamental model, decreased by 50 K (fundamental red edge, 
hereafter FRE) are reported in the third and fourth columns of Table 1, respectively. The 
adopted uncertainty of 50 K in the location of the IS boundaries stems from our adopted 



adopted in our previous theoretical papers on Cepheids ( 


Bono. Marconi. & Stelline 


wcrf 


1999; 


Fiorentino et al. 


2002; 


Marconi. Musella. & Fiorentino 


2005; 


Fiorentino et al. 


2007 



7 For models at 13 M Q we extrapolated the ML relation obtained for lower masses taking 
into account the evidence that, for higher metal content, models with 13 M follow the same 



ML relation as lower mass models. 
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e.g.) and RR Lyrae stars ( 


Bono et al. 


1997 




Vlarconi. Caputo. Di Criscienzo. & Castellani 


2003; 


Di Criscienzo. Marconi. & Caputo 


200^ 


)• 



For the first overtone mode, canonical pulsators show a stable oscillation only for 
stellar masses smaller than or equal to 7 M©, while we verified that the first overtone IS 
disappears increasing the mass by only 0.1 M Q . Noncanonical first overtone models present 
a stable limit cycle only for masses < 4.5M . The first overtone IS boundaries are reported 
in Table 2. 

Our theoretical results confirm that the longest period for first overtone pulsators 
(P FO max) increases as the metal content decreases, being of about 9 days at Z =0.0004 and 



of abo ut 6 days at Z=0.004 (IBono et al. 



2001 



20021 ). As already discussed by IBono et al. 



(120021 ). this behaviour is in marginal agreement with observational evidence in the Milky 
Way and in the Magellanic Clouds, and at variance with results from linear pulsation 
models that predict the opposite trend. 

The location in the HR diagram of fundamental and first overtone IS boundaries is 
shown in Figure 1, as obtained from canonical (middle panel) and noncanonical (upper 
panel) pulsation models, respectively. In this plot the solid lines represent the FBE 
and FRE, while the dashed lines are the boundaries of the first overtone instability 
strip. We notice that at fixed luminosity level, there is no significant difference in the 
effective temperature of the IS boundaries, as we move from a canonical to a noncanonical 
mass-to-luminosity (ML) relation. For canonical fundamental models, the bottom panel 
of Figure 1 shows the effect of a different assumption on the efficiency of convection, 
namely l/H p = 2.0 (dotted lines) instead of l/H p = 1.5. We notice that increasing the 
mixing-length parameter, the instability strip becomes na rrower, as already foun d in our 



previous investi gations of more metallic Cepheid models (iFiorentino et al 



RR Lyrae stars (IDi Criscienzo. Marconi, fc Caputo 



200J). 



2007J) and for 
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3.1. Comparison with the predictions at higher metallicities 



Figure 2 shows the comparison between the IS obtained in the present pa- 
per for fundamental canonical mode ls with l/H p = 1.5, and those obtained by 



Bono. Marconi. &: StellingwerJ (ll999f ) at metallicities traditionally assumed for the Magel- 
lanic Clouds and the Galactic Cepheids, namely Z = 0.004, Z = 0.008, and Z = 0.02. We 
notice that the location of the new boundaries confirms the trend already found at higher 



meta l content up to Z = 0.03 (see 



Fiorentino et al. 



2002 



Marconi. Musella. fc Fiorentino 



20051 ) that the predicted IS moves toward lower effective temperatures as the metallicity 
increases. In fact, the FBE for Z = 0.0004 is bluer than the one for Z = 0.004 because of 
the increased abundance of Hydrogen that produces an higher efficiency of the Hydrogen 
ionization region in driving pulsation. For the same reason, the FRE for Z = 0.0004 is 
cooler than the one for Z = 0.004 for masses lower than about 9M . At higher masses and 
luminosities (lower densities), the simultaneous effect of the lower metallicity, combined 
with the high efficiency of convection (that reduces the driving role of both Hydrogen and 
Helium ionization regions), prevails and produces a hotter red edge. 

However, the effect seems to be relatively small, so that no relevant implication on the 
zero point and on the slope of the theoretical PL relations are expected when moving from 
the typical metal contents of the Magellanic Clouds to the significantly lower abundance of 
IZwl8 (see below). 



4. An atlas of theoretical light curves for ultra-low metallicity Cepheids 



The bolometric light curves we have obtained for canonical fundamental pulsation 
models are shown in Figs. 3-8 for stellar masses ranging from 4 to 13 M Q . The light curve 
data for the complete set of models are available upon request from the authors. For stellar 
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masses ranging from 5 to 11 M the predicted light curves for l/H p = 2.0 are also shown 
for comparison. Confirmi ng our previous results on the effect of the convective efficiency on 



the pulsation properties ( IFiorentino et al 



2007 



Di Criscienzo. Marconi, fc Caputo 



2004J), 



the amplitude of the theoretical light curves decreases as l/H p increases from 1.5 to 2.0. 
An interesting feature shown by these models is the behaviour of the bump-Cepheids. As 
well known, Classical Cepheids in the period range 6 < P <16 days show a secondary 
maxi mum (bump) along both light and radial velocity curves. The Hertzsprung progress ion 



(HP; 



Hertzsprung 



1926 



Ledoux fc Walraven 



1958 



Bono. Marconi, fc Stellingwerl 



2000|) 



is the relationship linking the phase of the bump and the pulsation period. Inspection 
of Figure 5 suggests that at Z=0.0004 the transition from light curves with the bump 
on the decreasing branch to light curves with the bump on the increasing branch occurs 
at a pulsation period of ~ 13 days. This value is significantly longer th an the observed 



(IMoskalik. Buchler. fc Maromlll992l; IMoskalik et al. 



periods (IBono. Marconi, fc Stellingwerl 



2000 



Beaulieu 



19981 ) and predicted 



20001 ) at the center of the HP (P H p) for Galactic 



(« 10 d), LMC (« 10.5 d) and SMC (« 11 d) bump Cepheids. This theoretical 
result thus confirms the predicted an d observed trend of Php with metallicity (see also 



Marconi. Musella. fc Fiorentin 



mol 



20051 ). 



5. Fundamental relations for Z=0.0004 Cepheids 

The combination of the period- density relation for pulsating stars and the Stephan- 
Boltzmann law provides a strict relation between the pulsation period, the luminosity, the 
effective temperature and the mass. In Tables 3 and 4, we report the stellar masses, the 
luminosities, the effective temperatures, the periods and the predicted intensity-averaged 
mean magnitudes and colors (see below) for fundamental and first overtone models, 
respectively. A linear regression fit to the data reported in the first four columns of these 
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tables allows us to derive the following theoretical pulsation relations: 

log P = 10.754 + 0.921 log L - 0.770 log M - 3.336 log T e 

for fundamental models, and 

log P = 10.807 + 0.818 log L - 0.618 log M - 3.324 log T e 

for first overtone models, both showing an intrinsic dispersion of 0.01 mag. The bolometric 
light curves of the new models were transformed into the o bservational Johnson- Cousin s 



bands 



fll997al 



UBVRIJK) by means of the atmosphere models by ICastelli. Gratton. fc Kurucz 



bj) to derive the intensity-averaged mean magnitudes and colors reported in Tables 3 
and 4. Magnitude-averaged values were also computed and are available upon request from 
the authors. Static magnitude values have been derived and used to obtain the boundaries 
of the instability strip, at each chemical composition, in the various period-magnitude 
planes and, in turn, to construct synthetic multi band PL relations for fundamental 
canonical models (see Table 7). 

Linear regressions through the periods, magnitudes and colors reported in Table 3 
allow us to derive theoretical PLC and Wesenheit relations for fundamental Cepheids, both 
in the canonical case and in the noncanonical one. Since the two cases differ only by a 
luminosity offset (0.25 dex), the relations are the same in the two cases except for the value 
of the term logL/Lc, where L/Lc is the ratio between the luminosity and the canonical 
luminosity level for the corresponding chemical composition. The PLC has the form: 

< My > = a+/31ogP+7 (color) +5\ogL/L c , 

where a, ft, 7 and 5 have the values listed in Table 5 for each color. The Wesenheit relation 
has the form: 



< My > - 7 (color) = aw+Pw^ogP+dw^ogL/L 



c, 



13 



where aw, Pw, and Sw have the values listed in Table 6. 



The coefficients for the corresponding relations for first overtone pulsators are available 
upon request. Figures 9 and 10 show the comparison between models at Z=0.0004 (IZwl8) 
and Z=0.008 (LMC) in the PLC and Wesenheit planes, respectively, when the canonical 
ML relation is assumed. The magnitude differences at fixed period are mainly due to the 
effect of the different chemical composition on the predicted luminosity of each stellar 
mass. No relevant dependence on the assumed l/H p parameter is expected. The PLC does 
not depend on the instability strip topology and the Wesenheit relation has only a weak 
dependence, thanks to the inclusion of a color term with a coefficient given by the ratio 
between the t otal and the relative ex tinction and thus generally different from that of the 



PLC (see also 



Fiorentino et al 



2003). 



] 


following the same procedure adopted in our previous pai 


)ers ( 


Caputo, Marconi, & Musella 


2000 


Fiorentino et al. 


2002; 


Marconi. Musella. & Fiorentino 


2005 


), we have used the infor- 



mation on the boundaries of the IS and the pulsation r elation for fundamental mode 



Caputo. Marconi, fc Musella 



s to 



2000 



derive synthetic PL relations in the BVRIJK filters (see 
for details): 

Mj = A c + B c log P, 

where j is the photometric band. The coefficients A c and B c of the linear relations are 
reported in Tables 7 and 8 for canonical models with l/H p = 1.5 and 2.0, respectively, and 
in Table 9 for noncanonical models with l/H p = 1.5. (A nc and B nc ). Quadratic relations 
are also available, upon request from the authors. However, as discussed in the text, the 
nonlinear effect is very small at the chemical composition of IZwl8's Cepheids. 

Figure 11 shows the M^-logP distribution of canonical fundamental pulsators for the 
two different assumptions on l/H p (left and right panels) and for the noncanonical models 
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with I /Hp = 1.5 (middle panel). The solid lines represent the linear fits from Table 7. As 
already found for higher metal contents, the distribution of the pulsators becomes narrower 
when moving towards the infrared. We also find that the coefficients of the PL relations 
depend on both the ML relation and the value of the mixing length parameter (as also 
found from the comparison of Table 7 with Tables 8 and 9), but with an effect on the 
inferred distance lower than 10 %. 

The PL relations at Z=0.0004 are remarkably more linear than those obtained at 
higher Z with no clear evidence of a 



2008 



Sandage. Tammann. &: Reindl 



jreak at 10 days (see iNgeow. Kanbur fc Nanthakumar 



2009 



and references therein). In Figure 12 we show 
the comparison between the synthetic PL relations described here and those corresponding 
to the chemical compositions of Cepheids in the Magellanic Clouds (Z=0.004, Z=0.008 
Y=0.25), and in the Milky Way (Z=0.02 Y=0.28). The PL relations for Z=0.0004 are 
steeper than those at higher metallici ties, confirming the trend with metallicity predictec 



by our previous pulsation models (see 



Bono et al. 



2008 



Caputo. Marconi, fc Musella 



2000), 



but are very similar to those at Z=0.004 for periods longer than about 10 days. The 
difference in absolute magnitude between the predicted PL at Z=0.0004 and Z=0.004 is 
lower than 0.1 mag for periods between about 8 and 70 days. 



6. Theory versus observations 

In this section we compare the theoretical relations derived in the previ ous sections with 



the p roperties of the three Classical Cepheids we have identified in IZwl8 (IFiorentino et al. 



20101 ). It should be noted that the PL is a statistical relation, and it is not safe to apply 
it to a very small sample of Cepheids as that of IZwl8. For this reason, in the following 
we will focus only on the Wesenheit relations, that have much lower intrinsic scatter 
caused by the finite width of the instability strip, and on the light curve model fitting 
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technique. In particular, due to the unusually long periods of t wo of the 3 Cepheids 



in IZwl8 (namely, P=125 for V15, and 130.3 days for VI, see 



Fiorentino et al. 



2010 



bund 
for 



the identification of the variables), we will apply these two theoretic al methods only to th e 



bona fide classical Cepheid in IZwl8, namely star V6 (P=8.71 days, 



Fiorentino et al. 



2010h . 



In the two panels of Figure 13 we show the comparison between models and observations 
in the V, I Wesenheit versus period plane for two different assumptions on the ML relation. 
As these relations are reddening free, the shift required to superimpose the observed 
distributions to the theoretical relations directly gives the intrinsic distance modulus. By 
applying to V6 the canonical and noncanonical Period- Wesenheit relations in the V, I 
bands, we obtain two estimates of the distance modulus of the galaxy and the associated 
uncertainties, as labelled in each panel of Fig. 13. The two inferred values are consistent 
with each other within the errors and in a greem ent with our previous estimates of the 



distance to IZwl8 obtained in 



Aloisi et al 



(120071 ) . and in particular with the distance based 



on the TRGB (> Q = 31.30 ± 0.17 mag). 

An independent constraint on the distance to IZwl8 can also be obtained by applying 
the light curve model fitting technique to V6. This is a powerful tool to derive a direct 
estimate of the in trinsic stellar parameters, that has already been applied with success 



to both Galactic flNatale. M arconi 



1997 



Bono. Castellani. fc Marconi 



& Bono 



2002 



Keller fc Wood 



2002 



2008 1) and Magellanic f|Wood. Arnold, fc Sebo 



20061 ) Cepheids. Results 



from the application of the method to the observed V, I light curves of V6 are shown in 
Figs. 14 and 15, for isoperiodic sequences of models at fixed luminosity and at fixed mass, 
respectively. These models reproduce the observed curves within ± 0.13 mag simultaneously 
in V and I. In turn, we find that V6 is expected to have a stellar mass ranging from 
5.9M Q to 6.8M , log L/L & in the range from 3.64 to 3.69 dex and T e varying from w 
5800 K to « 6000 K. The corresponding average apparent distance moduli in V and I are 
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< pLv >=< fti >= 31.5 ± 0.2 mag, where the uncertainties take into account both the 
observational error and the intrinsic dispersion of the moduli obtained from the models 
shown in Figs. 14 and 15. By correcting the obtained apparent moduli for the extinction in 



the two bands ( adopting for IZwl 8 a 



consistent with 



Aloisi et al 



2007 and 



'oreground Galactic ex tinction E(B — V A )=0.032 mag, 



Fiorentino et al. 



20101), we find /x = 31.4 ± 0.2 mag, 



again in very good agreement with the distance modulus derived from the TRGB method 
Oo = 31.30 ±0.17 mag). 

As for the intrinsic stellar parameters inferred from the model fitting technique, we 
note that the mas s and luminosity estima ted for V6 are consistent with the evolutionary 



values obtained in 



Fiorentino et al. 



( 120101 ). According to the canonical ML relation adopted 
in this paper for a stellar luminosity log L/L Q = 3.68 ± 0.01 dex, we expect a stellar mass 
in the range from 6.0 to 6.5 M . Slightly smaller masses (down to 5.8 M Q ) can also been 
accepted, when invoking a moderate overshooting (5logL/L & <0.25 dex). On the other 
hand, stellar masses larger than ~ 6.5 M & seem to be excluded for the quoted luminosity 
level, at least on the basis of current stellar evolution models. Taking into account that 
the pulsation best fit model which minimizes the residuals between predicted and observed 
light curves both in V and in / bands, corresponds to a 6.5M with logL/L Q = 3.68 dex, 
that is very close to the canonical ML relation, we may conclude that the light-curve model 
fitting of V6 favours the canonical evolutionary scenario. 



7. Conclusions 

We updated our theoretical scenario for Classical Cepheids extending it to the ultra-low 
metallicity composition of the blue compact galaxy IZwl8, i.e. Z=0.0004. On the basis of 
these new models, we have derived and listed the coefficients for the key Cepheid relations 
(namely PLC, Wesenheit and synthetic PL relations) at Z=0.0004. 
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In addition, we find that: 



1. The predicted IS topology confirms our previous results on its dependence on the ML 
relation and assumed mixing-length parameter. 

2. The IS at Z=0.0004 is marginally bluer than that at Z=0.004. The small sensitivity 
to metallicity in this low Z regime leads to a leveling off of the metallicity effect on 
the synthetic PL relations. 

3. The predicted bolometric light curves show the already known dependence on 
effective temperature and mixing-length parameter. For bump Cepheid models 
the transition from light curves with the bump on the decreasing branch to 
light curves with the bump on the increasing branch occurs at a pulsation 



period of about 13 days. This value is significantly 



(IMoskalik. Buchler. &: Marom 



1992 



Moskalik et al. 



predicted periods (IBono. Marconi, fc Stellingwerf 



longer than both o 



Dserved 



2000 



Beaulieu 



19981 ) and 



20001) at the center of the HP (P HP ) 



for Galactic, LMC and SMC bump Cepheids. This o ccurrence confirms the predictec 



and o bserved trend of Pup with metallicity (see also 



2005|). 



Marconi. Musella. fc Fiorentino 



4. The synthetic PL relations for Z=0.0004 are more linear and steeper than those 
obtained at higher Z, but very similar to the Z=0.004 ones, at periods longer than 
about 10 days. The difference in absolute magnitude between the predicted PL at 
Z=0.0004 and Z=0.004 is lower than 0.1 mag for periods between about 8 and 70 
days. This result indicates a significant reduction of the metallicity effect on the 
Cepheid PL relations for Z < 0.004, thus supporting recent observational evidence of 
the universality of th e PL slope in galaxies with m etallicities significantly lower than 



in the LMC (see e.g. 



Pietrzynski et al 



2006 



20071 ). 



- 18 - 



5. The application to the IZwl8 Cepheid V6 of the Wesenheit relations derived with both 
canonical and noncanonical assumptions, provides two distance moduli consistent 
within the uncertainties with each other and with the distance to the galaxy inferred 
from the TRGB method. 

6. Also the model fitting of the light curves of V6 provides a distance modulus in good 
agreement with the TRGB distance. 

7. The intrinsic stellar properties obtained from the model fitting of V6 are consistent 
with the evolutionary estimates based on current sets of stellar models slightly 
favouring the canonical scenario. 

Finally we notice that, further observations of Cepheids in very metal poor galaxies 
are needed in order to clarify the nature of the very long period candidates discovered in 
IZwl8. If they were confirmed to be the extension of Classical Cepheids to higher masses, 
they would provide key constraints on pulsation models of massive stars. 

MM, IM, GC and MT acknowledge financial support from PRIN-INAF2005, PRIN- 
INAF2006 and COFIS ASI-INAF 1/016/07/0. Support for Francesca Annibali was also 
provided by NASA through grants associated with program GO- 10586. 
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Table 1: Predicted location of the IS boundaries for fundamental mode pulsators with 
Z=0.0004 and Y=0.24. 
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Table 2: Predicted location of the IS boundaries for first overtone pulsators with Z=0.0004 
and Y=0.24. 
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Table 3: Intrinsic stellar parameters, periods and intensity-averaged mean magnitudes and 
colors for fundamental models with Z=0.0004, Y=0.24, and l/H p = 1.5. 
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A portion of Table 3 is shown here for guidance regarding its form and content. The entire 
catalog is available in the electronic edition of the journal. 
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Table 4: As in Table 3, but for first overtone models at Z=0.0004, Y=0.24 and l/H p = 1.5. 
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Table 5: Coefficients of the multifilter PLC relations for fundamental mode models 
with Z=0.0004, Y=0.24, as a function of the adopted ML relation, in the form < 
My >=a+/31ogP+7(color)+<51ogL/L c 
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Table 6: Coefficients of the multi-filter Wesenheit relations for fundamental mode models 
with Z=0.0004, Y=0.24, as a function of the adopted ML relation, in the form < My > 
-'y(co\or)=a w +f3 w logP+5 w \og L/L c 
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Table 7: Coefficients of the linear synthetic PL relations based on canonical fundamental 
mode models with Z=0.0004, Y=0.24 and I /Hp = 1.5. Relations have the form Mj = 
A c + B c log P, where j is the photometric band. 
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B c 


a 


M B 


-0.76 


-2.81 


0.28 




±0.02 


±0.02 




My 


-1.10 


-3.04 


0.22 




±0.02 


±0.02 




Mr 


-1.36 


-3.14 


0.18 




±0.01 


±0.01 




Mj 


-1.67 


-3.24 


0.16 




±0.01 


±0.01 




Mj 


-2.03 


-3.37 


0.12 




±0.01 


±0.01 




Mr 


-2.36 


-3.50 


0.07 




±0.01 


±0.01 
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Table 8: Coefficients of the linear synthetic PL relations based on canonical fundamental 
mode models with Z=0.0004, Y=0.24 and l/H P = 2.0. Relations have the form Mj = 
A c + B c log P, where j is the photometric band. 



Magnitude 


A c 


B c 


a 


M B 


-0.97 


-2.80 


0.19 




±0.02 


±0.02 




My 


-1.27 


-3.03 


0.16 




±0.02 


±0.02 




Mr 


-1.51 


-3.13 


0.15 




±0.02 


±0.02 




Mj 


-1.79 


-3.23 


0.14 




±0.01 


±0.01 




Mj 


-2.12 


-3.36 


0.12 




±0.01 


±0.01 




M K 


-2.43 


-3.48 


0.10 




±0.01 


±0.01 
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Table 9: Coefficients of the linear synthetic PL relations based on noncanonical fundamental 
mode models with Z=0.0004, Y=0.24 and I /Hp = 1.5. Relations have the form Mj = 
A nc + B nc \ogP, where j is the photometric band. 



Magnitude 


A 


H nc 


a 


M B 


-1.01 


-2.54 


0.36 




±0.04 


±0.04 




My 


-1.22 


-2.85 


0.27 




±0.02 


±0.02 




Mr 


-1.42 


-2.98 


0.23 




±0.02 


±0.02 




Mj 


-1.68 


-3.10 


0.20 




±0.02 


±0.02 




Mj 


-1.98 


-3.25 


0.14 




±0.02 


±0.02 




Mr 


-2.26 


-3.40 


0.09 




±0.01 


±0.01 





3.9 



3.8 



3.7 



3.6 



log Te 



Fig. 1. — Location in the HR diagram of the fundamental (solid lines) and first overtone 
(dashed lines) IS boundaries for canonical (middle panel) and noncanonical (upper panel) 
models with l/H p = 1.5. In the lower panel the theoretical fundamental boundaries for 
canonical models with l/H p = 1.5 (solid lines) are compared with the same predictions at 
I /H p = 2.0 (dotted lines). 
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Fig. 2. — Comparison between the theroretical IS obtained in the present pa- 
per for fundamental canonical mo dels with l/H p = 1.5 and results obtained by 
Bono. Marconi, fc Stellingwerj fjl999f ) for Z = 0.004, Z = 0.008 and Z = 0.02. 
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Fig. 3. — Predicted bolometric light curves for canonical fundamental pulsation models with 
M = \M.. 
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Fig. 4. — Same as Figure 3, but for M = 5M . The dashed lines are the corresponding 
curves for l/H p = 2.0. 
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Fig. 5.- 



Same as Figure 4, but for M = 7M ( 
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Fig. 6.— Same as Figure 4, but for M = 9M ( 
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Fig. 7.— Same as Figure 4, but for M = 11M G 
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Fig. 8.— Same as Figure 3, but for M = 13M G 
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Fig. 9. — Comparison between theoretical PLC relations at 



Z=0.0004 and Z=0.008. 
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Fig. 10. — Comparison between theoretical Period- Wesenheit relations at Z=0.0004 and 
Z=0.008. 
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Fig. 11. — Synthetic multi-filter PL relations for canonical fundamental pulsators and two 
different assumptions on l/H p (left and right panels), and for noncanonical models with 
I /Hp = 1.5 (middle panel). 
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Fig. 12. — Comparison between the synthetic PL relations obtained in this paper, and the 
relations corresponding to the chemical composition of Cepheids in the Magellanic Clouds, 
and in the Milky Way. 
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Fig. 13. — Comparison between theoretical canonical (lower panel) and noncanonical (upper 



pan el) Period- Wesenheit r elations in the VI bands and data of the three Cepheids in IZwl8 



sec 



Fiorentino et al.l 



20101 ) . Distance moduli obtained from th e application of these relations 



to the bona fide Classical Cepheid V6 (jFiorentino et al. 



2010) are labelled. 



Fig. 14. — Model fitting of the light curves of V6 in the V and I band for a range of model 
masses and effective temperature (see figure's labels) at fixed luminosity level (see text). All 
the plotted models are able to reproduce the observed V,I light curves of V6 within 0.13 
mag (mean residual). - 44 - 
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Fig. 15. — Same as Fig. 14, but for models at fixed stellar mass. 
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